Dark spinor model with torsion and cosmology 
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We study cosmological consequences of the dark spinor model when torsion is included. Only some 
components of the torsion are allowed to be non- vanishing in homogeneous and isotropic cosmology, 
but there exist freedoms in the choice of these components which is consistent with the evolution 
equations. We exploit this and discuss several cases which can result in interesting cosmological 
consequences. Especially, we show that there exist exact cosmological solutions in which the Universe 
began its acceleration only recently and this solution is an attractor. This corresponds to a specific 
form of the torsion with a mild fine-tuning which can address the coincidence problem. 
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(N ; I. INTRODUCTION 

To search for the theoretical foundation of the current accelerating Universe [l[ has become one of the most 
fundamental problems in modern cosmology. An unknown source of energy which is repulsive in nature referred to as 
$H . the dark energy is one of the most promising candidates and various proposals for its origin have been put forward Q . 
But so far, relatively little attention has been paid to the possibility of fermionic sources in cosmology even though 
it has been known for some time that spinor fields could play some important roles in the evolution of the Universe 
CsJ ■ For example, in Ref. [|| it is shown that a spinor field can accommodate any desired behavior of its energy 

density if an appropriate self-interaction of the spinor field is introduced and it is suggested that it can account for 
the dark energy and allows for realistic cyclic non-singular solutions. In Ref. Q it is investigated whether fermionic 
sources could be responsible for the accelerating periods during the evolution of a universe after a matter field would 
guide for the decelerating period. In Ref. [6| the authors have shown that it is possible to simulate perfect fluid and 
dark energy by means of a nonlinear spinor field. 

Recently, the so-called ELKO spinor, which is a new spin-half field defined as the eigenspinor of the charge conju- 
gation operator with mass dimension one, was introduced in Ref. [7[. They have many special properties, one being 
that their dominant coupling is via the gravitational field and therefore are naturally dark. It has been proposed to 
be a candidate of dark matter Q and dark energy H and many cosmological aspects of the ELKO spinor model have 
been investigated M. In addition, dark spinor was shown to be a potential candidate to probe non-trivial topolo- 
gies in spacetime [10j | and to be incorporated in a natural extension of the Standard Model by representing a mass 
^ ■ dimension-transmuting operator between Dirac and ELKO spinor [llj . 

In Ref. [IH, however, it was pointed out that the construction of ELKO spinors itself implicitly violates Lorentz 
invariance. The authors have proposed a non-local but Lorentz invariant version of the dark spinors. Also it is pointed 
out that some crucial errors of neglecting the spin connection part of the fermionic sector in the computations of stress 
energy tensors were made in the previous literatures on ELKO models. The correct expression for the ELKO spinor 
field was applied to dark energy and the result showed the existence of the de Sitter acceleration of the Universe. 
Another interesting consequences of dark spinor in cosmological aspects is that the equation of the state ui for the dark 
spinor can be negative and can cross the phantom-divide even without any phantom-like matter carrying negative 
kinetic energy [13|. It converges to a parameter value within the bound produced by using WMAP data hj as was 



discussed in Ref. 
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In spite of these nice features, further dynamical analysis of the ELKO model following this work using the correct 
form of the stress energy momentum tensor revealed that scaling attractor solutions do not exist in this model [l5[ , 
and it has the difficulty of addressing the coincidence problem. This aspect of the dark spinor model draws our 
attention and motivates to extend to the case in which torsion (l6j is included in the theory. First, recall that when 
the cosmological principle is extended to the general relativity theory with torsion, only the time component of the 
trace and totally anti-symmetric components of the torsion tensor can be nonzero [l7j . For ordinary Dirac spinor, 
it is well-known that it interacts only with the totally anti-symmetric components of the torsion tensor and the spin 
connection part does not contribute to the total energy density (l8| . On the other hand, the dark spinor interacts 
with all components of the torsion and we investigate the cosmological consequence of the torsion in the dark spinor 
model. In particular, we find the time component of the non- vanishing trace vector could provide a novel feature of 
explaining the coincidence problem. That is, with a mild fine-tuning of the parameters, we show that there exist a 
cosmological solution in which the Universe began its acceleration only recently and that the solution is an attractor. 

The paper is organized as follows: In Sec. II, we give a brief summary of the spinor field in curved background with 
torsion in order to set the notation. In Sec. Ill, the action of the dark spinor is given and the equations of motions 
are calculated. In Sec. IV, we discuss various cosmological solutions focusing on the solution which can address the 
coincidence problem and stability analysis of this solution is performed. Sec. VI contains conclusion and discussions. 

II. SPINOR IN CURVED BACKGROUND 

In this section, we give a brief summary of the spinor field in curved background which is essential for our analysis. 
More detailed descriptions can be found in Ref. [l9(. The gamma matrices are given by 



with {7°, 7 b } = -2r/ ab , r] ab = diag.(-l, 1, 1, 1) 
We introduce the metric and tetrad via 



9^ = e"4??a6- ( 2 ) 



Imposing the metricity V^i?^ = 0, the connection is given by 

r£„ = {^}-V (3) 

where {} is the Christoffel symbol and the contortion K l±v p is given by the torsion tensor S IIV P = Ff , via [l|| 

Kfj,v P — ~ ip^v P + S P fiv + S P vn) ■ (4) 
The covariant derivative of the spinor and its dual are given by 

= - r M v (5) 

and 

V^ = c^ + #V, ( 6 ) 

where T M is the connection on the spinor. 

We introduce the spin connection by imposing the metricity for the tetrad 

o = v (1 e» = a (1 e:-r^;+< 6 et, (7) 

where the spin connection satisfy lj a b = — oj h ° and 

^ ab = e a v (d^ + T^ p e b "). (8) 
We also impose the covariant constancy of the gamma matrices; 

= V, t7 p - d,Y + T^-f - [IV, 71- (9) 
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From the eqs Q and ([9]), we obtain 

IV = J< b 7a7b + cl. (10) 

We consider only the case when the constant c is equal to zero. 

We can decompose the contortion tensor ^ into a traceless part and trace [20l ]: 

K^^K^-^SPSu-g^S"), (11) 

where K^ p is the traceless part, K ^ = and S„ is the trace of the torsion tensor, S v — . Making use of (|lll) . 
we can write curvature scalar as follows: 

R = R- 4V M S" - - k vpa K av P. (12) 

In the right hand side of the above equation, scalar curvature and the covariant derivative are calculated with respect 
to the Christoffel symbol. 



III. TORSION DARK SPINOR MODEL 



Let us consider the dark spinor action of the form 



S = d i xy / —g 



where R is the scalar curvature including the torsion piece and V is the potential for the dark spinor (ELKO) field ip 
and its dual ip. S m is the matter field action. 

To discuss the Friedman cosmology in the fiat Robertson- Walker space-time, consider a metric of the form 



ds 2 = -N 2 dt 2 + a 2 (i)cfeW, 
where a(t) is the scale factor of our three dimensional universe. We assume 



(13) 



(14) 



with a homogeneous real scalar field (j)(t) and a constant spinor £ such that = 1 but V M £ ^ 0. It is to be pointed 
out that that the spinor nature of the dark spinor is remnantal in the real scalar field 4>{t) defined in ([M)) and the 
ensuing evolution equations contain the contributions from the spin connection which are absent in the case of a 
genuine scalar field. It was shown that these new contributions supply fundamentally different aspects in the dark 
spinor cosmology 12 1. 

For isotropy and homogeneity of Universe, one assumes the following non- vanishing components of the torsion (l7j 



So 



5'. 



.in 



h(t)/2, V = 5 b .; ] = f(t)/3ej 



(15) 



Eqs. (jl3|) and (fT5j) give the following connection components: 

„n N „ n ab 



N 



r° 



N 



-t) T l 



Toi — Fjk — f e ijk, 

a 



with b = a + ah. Using the above components, we obtain the scalar curvature 



R=Q[ 



1 d 



b 

aN 



Note that the above result ([TT]) agrees with Eq. ^21 up to a total derivative term. From the metric (fl 
is given by 



(N,a,a,a), 



1111 

TV' a' a' a 



(16) 

(17) 
the tetrad 

(18) 
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Using (|16p , the only non- vanishing components of the spin connection ((5J are given by 



3 ~ u *° - 



Then, we have from eq.(|T 



x* ""a; 1 3 



T n = 0, T„i = 



b f 



Putting everything together, we have (k = 1) 



S 



dt 



1 

N 



-3ad 2 + 3a 3 h 2 + -a 3 ^ 2 



-ah 1 



N 3a f + -af 2 <t> 2 + a 3 V 



The equation related with the energy density coming from the variation with respect to N is given by 



3H 2 = -<j> 2 + V + -H 2 <j) 2 + -Hh cf> 2 + 3 i I 
2 8 4 



1 



3 1 



1 



11 + 

n ^ 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



Note that the above equations (|2"2"|) . ([M]). ([2^)1 arc different from the equations one would obtain by directly using 
Gfj,i, = TJJl, where G M „ is the Einstein tensor calculated with the connections including the torsion as in Eq. (fT()]) 
pH [22l | . The discrepancy comes from the fact that when torsion is included, the variation SR^ of the Ricci tensor in 
calculating the equations of motion from the scalar curvature is not a total derivative and cannot be neglected as in 
the case without torsion. Therefore, G p „ — T™ is not justified in the case when torsion is included. Differentiation 
Eq. (|2"2"|) with respect to time and using Eqs. (|24l) and (|2"3"1) . we obtain the following consistency equation: 



and the one related with the pressure coming from the variation with respect to a is 



1 



1 d 



2H - 3H 2 = -4> 2 - V ~ -H 2 4> 2 -- A j t [{H + h)4> 2 ]+i[l + -4> 2 )h 2 + i[l + -<j ) 2 )^+p. 



2^ 8 
The scalar field equation is given by 



= 4> + 3H<j) + V'{4>) --i{H + h) 2 - f — 



1 



2h + QHh +-(H + h), 



1 



= 3/i 

which serves as a dynamical equations for the torsion. 



2 + -H(H + h) <j> 2 



IV. COSMOLOGICAL SOLUTIONS 



In this section, we discuss diverse cosmological solutions depending on the torsion. 



A. Torsion as a dark matter candidate 



We first analyze the pure torsion case where the dark spinor field is absent, (f> = 0. Let us also assume that 
Pm = Pm = 0. If we define ph = 3h 2 and Pf = 3/ 2 /a 2 , we have from Eqs. (|22|) and (|23| 

3H 2 = p h + p f , -2H-3H 2 = Ph -l Pf . (26) 

The second equation says that ph = Ph, whereas pf = —pf/3. Eq. ([25]) gives the continuity equation for the torsion 
energy density: 

p h + 6Hp h + p f +2Hp f = 0. (27) 
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Let us suppose that ph and pf are separately conserved. Then, ph oc 1/a 6 corresponding to the equation of state 
uih = 1 an d / = fo is a constant corresponding to u)f — —1/3 as expected. However, if we assume Pf = rp^ with a 
constant r > and p t = pn + pf = (1 + and the pressure p t = ph — Pf/3 = (1 — r/3)ph, Eq. ff2T|) becomes 

/> 4 + 3Aff A = 0, (28) 
3 + 3r 

Therefore, p t can describe an ideal fluid with the range of the equation of state given by —1/3 < uj t = < 1. 

When r — 3, it describes the cold dark matter. It is to be commented that the relation pf = rp^ is very ad hoc, 
nevertheless nothing forces separate conservation as far as continuity equation (|2"T|) is concerned. 



B. Torsion and cosmological constant 



Next, wc consider the cases / = 0. Let us define a quantity 

X = 2h+^{H + h)(j> 2 . (29) 

We assume that there is no interaction between the matter and spinor dark energy with torsion, and the matter and 
the dark spinor are separately conserved. 

p m + 3H(p m + Pm ) = 0. (30) 

Then from Eq. (|25|) . we have 

3h (x + ZHX^j = 0. (31) 
We consider h ^ and then, Eq. (13T1) can be immediately integrated to give 

or 

Eq. ([23]) by using (J22J) becomes 

- 2 f H + h] - 4> 2 + 3(H + h)X + p m + Pm . (32) 



Also, Eqs. (|22l) and (IM1) can be rewritten as 



3H 2 = \j> 2 + V + ^ (X - 2hf + 3h 2 + p m , (33) 



O = + 3ff0 + V'(0)-12 (X J k)2 , (>=4l) (34) 



and 



We have four equations (|3Tj) . (I32|) . ([33|) . and (|34[) . These equations are not independent, and we have three equations 
to be solved including Eq. (f53"|). 

The torsion component h which was introduced for homogeneity and isotropy via (1151) is an external input as long as 
it satisfies the evolution equations. We exploit this arbitrariness to assume an existence of the functional dependence 
of h(t) = h{4>) and discuss cases where some analytic properties of the above equations can be investigated. 

Let us discuss the case where p m — p m = 0. Under the further simplifying assumption of c\ = 0, some exact 
cosmological solutions can be obtained in a couple of cases which describe the accelerating Universe. First, X = 
and we have 



H=-h[l + — \. (35) 



We are interested in the case of ft. < 0. From Eq. (|32[) . we have 



16 

= -~(ht<t>-2h), (36) 
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which yields from Eq. (|33l) the potential of the form 

y=**£(8 + ^)-^(*,^-2&) a . (37) 

One can check that Eq. flMJ) is satisfied with Eqs. (|3"5]l. (j31)l) and ([37]) . 

Let us discuss a couple of cases where exact cosmological solution can be obtained. The first case is when the 
potential is of the form 

V{(f>) = ^toV + K). (38) 

It corresponds to the choice of h oc cf> 2 . Comparing with the potential (|3"T|) . we find they coincide if we choose 
h = ±rruj) 2 /-\/48- Moreover, the cosmological constant Vb is given by 

V = 4to 2 . (39) 



One can check that Eq. (|34l) is also satisfied. From Eq. (f3T>)) , we have = 0, and <f> — 4>q. Choosing h — —m<p 2 /v48, 
the Hubble parameter is given by 

tf = ^(8 + ^)> (40) 



for the de Sitter accelerating solution. A couple of comments are in order: First, <j) need not be the minimum of 
the potential (|38[) for its stable equilibrium point and can stay away from zero. This can be seen dynamically from 
Eq. ([Ml) as follows: The restoring force coming from V' is exactly canceled by the repulsive force coming from the 
torsion if h cx cj) 2 . Therefore, the scalar field undergoes a damped motion without any force acting, and this solution 
corresponds to the solution in which if it takes initial value </>o, it stays there forever. Also, we note the difference with 
an ordinary quintessence model with mass term and cosmological constant. In this case, the <fi = is dynamically 
preferred, and the acceleration can be driven by the cosmological constant. However, the crucial difference is that 
there is a priori no relation between the cosmological constant and the mass of the scalar field, whereas in the dark 
spinor model with torsion, there is a relation given by eq. (1391) . This implies that the smallness of the cosmological 
constant is directly interwoven with the ultra light scalar field. 

C. Torsion and coincidence problem 

Another case of interest is when h is given by 

h = -|c| cxp(-\c/)/2)<j> 2 (X > 0). (41) 
In this case, we have potential of the form 

V = Ke" A * [cj) 2 + c] , V, = 24c 2 , = 8 - ^A 2 . (42) 



4> = A*e- A< ^ 2 (A* =8|c|A), (43) 



Eq. (J36J) gives 



which can be integrated to yield 

cj)(t) = hnt + A. (44) 
A 

We choose the integration constant A of the above equation to be zero which fixes a relation 4|c|A 2 = 1. The Hubble 
parameter H by using Eq. (|35l) and its first time derivative arc given by 



1 



(lnt) ; 
~2X r 



XH 2 



1 - 



kit (hit) 



X 2 2X 2 



(45) 
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To calculate the scale factor, let us define — d\na/d(j>. Then, Eq. (|35|) can be readily integrated to give 

a((f>) = a^ei^^ = 



To study the existence of the accelerating solution, we first check 



H + H 



1 

t 2 



A 2 



1 



2 

A 2 



( lnt 
' A 4 " 



1 



A\ (lnt) 2 (lnt) 4 
A 2 J + 



(46) 



(47) 



We see that at early times when A 3 / 2 < lnt < A 2 (A > 1), Eq. (|4"7} describes a decelerating universe. This is the 
epoch when the first term in H is dominant. When lnt > A 2 the (lnt) term begins to become dominant gradually 
and at later times , H 2 term becomes dominant. The acceleration would begin around lnt QCC ~ A 2 . If we adjust the 
parameter A to be ~ 12 we get the late time acceleration occurring around t acc ~ 10 61 t p i, which roughly corresponds 
to the current age of the Universe. This could explain the coincidence problem. The dependence of t acc on the 
numerical value of A is depicted in Fig. 1. Comparing with the A — CDM which requires a fine-tuning of the order 
of ~ 1CP 120 for the cosmological constant, it is of 0(1) and poses no fine-tuning issue. 



400 P 




~ 12, t acc ~ I0 m s) 



10 

A 



15 



20 



FIG. 1: The relation between t acc and the parameter A in the Planck unit. The current acceleration can be achieved with 
A « 12. 



We perform a stability analysis of the solution (|44l) and (1461) with A = 0. Variation of the Eq. (|34|) with X = 

gives 



= 50 + 3H50 + 3<j>6H + V"(<p)S<p + 144^^ - 96-^tit 



h 2 



SH can be eliminated through Eq. (|55|) via 

5H 

Using Eqs. (|gT). p2|) . and we find that Eq. g8]) reduces to 

3 



-/i'(l + — ) + 16— 



+ A (2A 2 + (lnt) 2 ) S4>+^ (/(A) + (lnt) 2 ) S<j> = 0, 



with /(A) being 



/(A) = — -A 4 + 2A 2 — -. 
J\ ) 4 3 



(48) 



(49) 



(50) 



(51) 



A numerical analysis of Eq. ([50)) shows the attractor behavior of the solution, as is shown in Fig. 2. 
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FIG. 2: Attractor behavior: The initial condition is given at Ni = Inti ~ 130 which corresponds approximately to the 
decoupling time. The red curve is with 4>(Ni) = 0.005, <j>{Ni) = 0.01, blue with 4>{Ni) = -0.001 ,4>(Ni) = 0.01, green with 
tj>(Ni) = -0.005 , <j>(Ni) = -0.01, and pink : <j>{Ni) = 0.001 , (j>(Ni) = -0.01. 

A dynamical estimate how the solution becomes an attractor can be given. To check the stability in the asymptotic 
region where lni >> A 2 , <f> >> 2 A, we first note that from Eq. (|42p 



V"56 ~ 24c z e 
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xn- 



i\nt) 2 8<t>, 



(52) 



and among the last three terms of Eq. (|4"8")l . V"(<j>) term dominates. This is a positive quantity and provides a 
restoring force in the asymptotic region. Also, from Eq. (|49)) . we have 8H ~ —h'5(f>, which yields using Eq. (|4T|) 

(53) 



35Hj> ~ = -^(lnt) 2 ^. 

This term acts as a source of a repulsive force, but their magnitude is smaller then the attractive force of Eq. ([5 
Over all, in the asymptotic region, the perturbation equation becomes 



3 (lnf) 2 • 3 {\nt)< 



A 4 t A 4 t 2 

It can be easily checked that the perturbation decays as 

(lnf) Q 



0. 



t ' 



(54) 



(55) 



and the solution is an attractor in the asymptotic region. 



V. CONCLUSION 



We have examined the possible cosmological role that the torsion could play in the dark spinor model. Freedom 
exists in the choice of the torsion, which is consistent with the continuity equation. This has been exploited for the 
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choice of diverse forms of the torsion, which resulted in the various cosmological consequences. In the pure torsion 
case without dark spinor, torsion could play the role of the dark matter. When the torsion is proportional to (j> 2 , we 
have an accelerating de Sitter Universe. In particular, we found an exact cosmological solution in which the Universe 
began its acceleration only recently and it is an attractor. This corresponds to a specific form of the torsion with a 
mild fine-tuning and it can address the coincidence problem. 

In this work, the torsion was assumed to be a kind of external matter field whose dynamics is dictated by the con- 
tinuity equation. Even though this does not confront with any theoretical inconsistency, their choices were restrictive 
to certain forms and could be regarded as rather ad hoc. In the more fundamental approach like the Einstein-Cartan- 
Sciama-Kibble theory of gravity [l6| where the torsion tensor is regarded as a dynamical variable from the beginning, 
the equations of motion coming from the variation with respect to the torsion generates a spin-spin interaction equa- 
tion. In the Dirac case, this interaction is significant only at extremely high densities and it was shown that such 
an interaction averts the unphysical big-bang singularity, replacing it with a cusp-like bounce at a finite minimum 
scale factor (23[. This approach could be extended to the ELKO case [24|. In Ref. (25[, the most general ELKO 
matter fields including the dynamical torsion was constructed, whose cosmology leads to anisotropic universes. To 
the knowledge of the authors, the relation between the two approaches has not been explored, and it seems that at 
present only the former approach could provide the novel cosmological feature of addressing the coincidence problem 
for the dark spinor. 

One interesting feature of the potential (|4"2"j) is that in contrast with the pure exponential potential in the quintessence 
model where the parameter A is restricted to be within a certain range to produce the attractor solution, the solution 
with the torsion is an attractor for each value of A > 0. Moreover, the term e~ x ^/ 2 4> 2 is reminiscent of the Albrccht- 
Skordis potential (2(| which appears in the low energy limit of M theory. The presence of the torsion is a key ingredient 
and possible implications of this aspect deserve further investigations. 

We conclude with the following remarks. As shown in Eq. (159")) of Sec. IV B, there is a relation between the 
cosmological constant and the mass of the cosmological particle cj) in our dark spinor model with torsion, which might 
represent the main component of our present universe and is about 10~ 63 <?, as in the massive vector model [13] ■ When 
the analysis is extended to the case including a conformal coupling of dark spinor fields to gravity, it will provide 
with a richer variety of cosmological consequences than the case without torsion [28|. Last but not least, possible 
implications of the experimental limits on torsion [29( in the dark spinor case need to be explored. 
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